Search for: All records

The Faraday rotation measure (RM) is a commonly used tool to trace electron number density and magnetic fields in hot accretion flows, particularly in low-luminosity accreting supermassive black holes. We focus on the nuclear region of M87, which was observed at 230 GHz (1.3 mm) by the Event Horizon Telescope in 2019. It remains unclear whether this emission originates from the accretion flow, the jet base, or both. To probe the presence of an accretion flow, we explore the scenario where the linearly polarized emission from the counter jet, visible at 43 GHz (7 mm), is Faraday-rotated by the accretion flow. We calculate theoretical predictions for counter-jet polarization using analytical and numerical models. In all cases, we find a Faraday-thick flow at 43 GHz (7 mm), with RM ∼ 10⁶rad m⁻², and a polarization angle that follows a linear relationship with wavelength squared, consistent with external Faraday rotation. The more realistic model, which includes turbulence and magnetic field fluctuations, predicts that the polarization pattern should be time-dependent, and that the counter-jet emission is depolarized due to Faraday depth fluctuations across the accretion flow. Despite the Faraday thick regime and strong depolarization, the linear relationship persists, enabling us to constrain the flow’s physical properties. Comparing the counter-jet and forward-jet linear polarization states should enable detection of M87’s accretion flow and provide lower limits on electron density, magnetic field strength, and mass accretion rate. 

Abstract In certain scenarios, the accreted angular momentum of plasma onto a black hole could be low; however, how the accretion dynamics depend on the angular momentum content of the plasma is still not fully understood. We present three-dimensional, general relativistic magnetohydrodynamic simulations of low angular momentum accretion flows around rapidly spinning black holes (with spina = +0.9). The initial condition is a Fishbone–Moncrief (FM) torus threaded by a large amount of poloidal magnetic flux, where the angular velocity is a fractionfof the standard value. Forf= 0, the accretion flow becomes magnetically arrested and launches relativistic jets but only for a very short duration. After that, free-falling plasma breaks through the magnetic barrier, loading the jet with mass and destroying the jet–disk structure. Meanwhile, magnetic flux is lost via giant, asymmetrical magnetic bubbles that float away from the black hole. The accretion then exits the magnetically arrested state. Forf= 0.1, the dimensionless magnetic flux threading the black hole oscillates quasiperiodically. The jet–disk structure shows concurrent revival and destruction while the gas outflow efficiency at the event horizon changes accordingly. Forf≥ 0.3, we find that the dynamical behavior of the system starts to approach that of a standard accreting FM torus. Our results thus suggest that the accreted angular momentum is an important parameter that governs the maintenance of a magnetically arrested flow and launching of relativistic jets around black holes. 

Abstract The presence of a strong, large-scale magnetic field in an accretion flow leads to the extraction of the rotational energy of the black hole (BH) through the Blandford–Znajek (BZ) process, believed to power relativistic jets in various astrophysical sources. We study rotational energy extraction from a BH surrounded by a highly magnetized thin disk by performing a set of 3D global GRMHD simulations. We find that the saturated flux threading the BH has a weaker dependence on BH spin, compared to highly magnetized hot (geometrically thick) accretion flows. Also, we find that only a fraction (10%–70%) of the extracted BZ power is channeled into the jet, depending on the spin parameter. The remaining energy is potentially used to launch winds or contribute to the radiative output of the disk or corona. Our simulations reveal that the presence of a strong magnetic field enhances the radiative efficiency of the disk, making it more luminous than its weakly magnetized counterpart or the standard disk model. We attribute this excess luminosity primarily to the enhanced magnetic dissipation in the intra-ISCO region. Our findings have implications for understanding X-ray corona formation and BH spin measurements, and interpreting BH transient phenomena. 

Abstract We carry out idealized three-dimensional general-relativistic magnetohydrodynamic simulations of prograde, weakly magnetized, and geometrically thick accretion flows where the gas distribution is misaligned from the black hole (BH) spin axis. We evolve the disk for three BH spins:a= 0.5, 0.75, and 0.9375, and we contrast them with a standard aligned disk simulation witha= 0.9375. The tilted disks achieve a warped and twisted steady-state structure, with the outer disk misaligning further away from the BH and surpassing the initial 24° misalignment. However, closer to the BH, there is evidence of partial alignment, as the inclination angle decreases with radius in this regime. Standing shocks also emerged in proximity to the BH, roughly at ∼6 gravitational radii. We show that these shocks act to partially align the inner disk with the BH spin. The rate of alignment increases with increasing BH spin magnitude, but in all cases is insufficient to fully align the gas before it accretes. Additionally, we present a toy model of orbit crowding that can predict the location of the shocks in moderate-to-fast rotating BHs, illustrating a potential physical origin for the behavior seen in simulations—with possible applications in determining the positions of shocks in real misaligned astrophysical systems. 

Abstract Magnetic reconnection is ubiquitous in astrophysical systems, and in many such systems the plasma suffers from significant cooling due to synchrotron radiation. We study relativistic magnetic reconnection in the presence of strong synchrotron cooling, where the ambient magnetization,σ, is high and the magnetic compactness,ℓ_B, of the system is of order unity. In this regime,e^±pair production from synchrotron photons is inevitable, and this process can regulate the magnetizationσsurrounding the current sheet. We investigate this self-regulation analytically and find a self-consistent steady state for a given magnetic compactness of the system and initial magnetization. This result helps estimate the self-consistent upstream magnetization in systems where plasma density is poorly constrained, and can be useful for a variety of astrophysical systems. As illustrative examples, we apply it to study the properties of reconnecting current sheets near the supermassive black hole of M87, as well as the equatorial current sheet outside the light cylinder of the Crab pulsar. 

Abstract Quasars show a remarkable degree of atomic emission-line broadening, an observational feature which, in conjunction with a radial distance estimate for this emission from the nucleus, is often used to infer the mass of the central supermassive black hole. The radius estimate depends on the structure and kinematics of this so-called broad-line region, which is often modeled as a set of discrete emitting clouds. Here, we test an alternative kinematic disk-wind model of optically thick line emission originating from a geometrically thin accretion disk under Keplerian rotation around a supermassive black hole. We use this model to calculate broad emission-line profiles and interferometric phases to compare to GRAVITY data and previously published cloud modeling results. While we show that such a model can provide a statistically satisfactory fit to GRAVITY data for quasar 3C 273, we disfavor it as it requires 3C 273 be observed at high inclination, which observations of the radio jet orientation do not support. 

Abstract We explore reprocessing models for a sample of 17 hypervariable quasars, taken from the Sloan Digital Sky Survey Reverberation Mapping project, which all show coordinated optical luminosity hypervariability with amplitudes of factors ≳2 between 2014 and 2020. We develop and apply reprocessing models for quasar light curves in simple geometries that are likely to be representative of quasar inner environments. In addition to the commonly investigated thin-disk model, we include the thick-disk and hemisphere geometries. The thick-disk geometry could, for instance, represent a magnetically elevated disk, whereas the hemisphere model can be interpreted as a first-order approximation for any optically thick out-of-plane material caused by outflows/winds, warped/tilted disks, and so on. Of the 17 quasars in our sample, 11 are best-fitted by a hemisphere geometry, five are classified as thick disks, and both models fail for just one object. We highlight the successes and shortcomings of our thermal reprocessing models in case studies of four quasars that are representative of the sample. While reprocessing is unlikely to explain all of the variability that we observe in quasars, we present our classification scheme as a starting point for revealing the likely geometries of reprocessing for quasars in our sample and hypervariable quasars in general. 

ABSTRACT Sgr A* exhibits flares in the near-infrared and X-ray bands, with the luminosity in these bands increasing by factors of 10–100 for ≈60 min. One of the models proposed to explain these flares is synchrotron emission of non-thermal particles accelerated by magnetic reconnection events in the accretion flow. We use the results from particle-in-cell simulations of magnetic reconnection to post-process 3D two-temperature GRMHD simulations of a magnetically arrested disc (MAD). We identify current sheets, retrieve their properties, estimate their potential to accelerate non-thermal particles, and compute the expected non-thermal synchrotron emission. We find that the flux eruptions of MADs can provide suitable conditions for accelerating non-thermal particles to energies γe ≲ 106 and producing simultaneous X-ray and near-infrared flares. For a suitable choice of current-sheet parameters and a simplified synchrotron cooling prescription, the model can simultaneously reproduce the quiescent and flaring X-ray luminosities as well as the X-ray spectral shape. While the near-infrared flares are mainly due to an increase in the temperature near the black hole during the MAD flux eruptions, the X-ray emission comes from narrow current sheets bordering highly magnetized, low-density regions near the black hole, and equatorial current sheets where the flux on the black hole reconnects. As a result, not all infrared flares are accompanied by X-ray ones. The non-thermal flaring emission can extend to very hard (≲ 100 keV) X-ray energies. 

ABSTRACT Magnetically arrested accretion discs (MADs) around black holes (BHs) have the potential to stimulate the production of powerful jets and account for recent ultra-high-resolution observations of BH environments. Their main properties are usually attributed to the accumulation of dynamically significant net magnetic (vertical) flux throughout the arrested region, which is then regulated by interchange instabilities. Here, we propose instead that it is mainly a dynamically important toroidal field – the result of dynamo action triggered by the significant but still relatively weak vertical field – that defines and regulates the properties of MADs. We suggest that rapid convection-like instabilities, involving interchange of toroidal flux tubes and operating concurrently with the magnetorotational instability (MRI), can regulate the structure of the disc and the escape of net flux. We generalize the convective stability criteria and disc structure equations to include the effects of a strong toroidal field and show that convective flows could be driven towards two distinct marginally stable states, one of which we associate with MADs. We confirm the plausibility of our theoretical model by comparing its quantitative predictions to simulations of both MAD and SANE (standard and normal evolution; strongly magnetized but not ‘arrested’) discs, and suggest a set of criteria that could help to distinguish MADs from other accretion states. Contrary to previous claims in the literature, we argue that MRI is not suppressed in MADs and is probably responsible for the existence of the strong toroidal field. 

Abstract We use the public code ebhlight to carry out 3D radiative general relativistic magnetohydrodynamics (GRMHD) simulations of accretion on to the supermassive black hole in M87. The simulations self-consistently evolve a frequency-dependent Monte Carlo description of the radiation field produced by the accretion flow. We explore two limits of accumulated magnetic flux at the black hole (SANE and MAD), each coupled to several subgrid prescriptions for electron heating that are motivated by models of turbulence and magnetic reconnection. We present convergence studies for the radiation field and study its properties. We find that the near-horizon photon energy density is an order of magnitude higher than is predicted by simple isotropic estimates from the observed luminosity. The radially dependent photon momentum distribution is anisotropic and can be modeled by a set of point-sources near the equatorial plane. We draw properties of the radiation and magnetic field from the simulation and feed them into an analytic model of gap acceleration to estimate the very high energy (VHE) γ-ray luminosity from the magnetized jet funnel, assuming that a gap is able to form. We find luminosities of $$\rm \sim 10^{41} \, erg \, s^{-1}$$ for MAD models and $$\rm \sim 2\times 10^{40} \, erg \, s^{-1}$$ for SANE models, which are comparable to measurements of M87’s VHE flares. The time-dependence seen in our calculations is insufficient to explain the flaring behaviour. Our results provide a step towards bridging theoretical models of near-horizon properties seen in black hole images with the VHE activity of M87.